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C^ , Abstract 

We consider a nonrelativistic quantum charged particle moving on a 
plane under the influence of a uniform magnetic field and driven by a 
periodically time-dependent Aharonov-Bohm flux. We observe an ac- 
celeration effect in the case when the Aharonov-Bohm flux depends on 
\G . time as a sinusoidal function whose frequency is in resonance with the 

cyclotron frequency. In particular, the energy of the particle increases 
linearly for large times. An explicit formula for the acceleration rate 



^— N . is derived with the aid of the quantum averaging method, and then it 

is checked against a numerical solution with a very good agreement. 

Keywords: electron-cyclotron resonance, Aharonov-Bohm flux, quan- 
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1 Introduction 

The problem of acceleration in physical systems driven by time-periodic ex- 
ternal forces, both in classical and quantum mechanics, has a rather long 
history though the results in the latter case are much less complete. One of 
the most prominent examples which initiated a lot of efforts in this field is 
the so called Fermi accelerator. On the basis of a theory due to Fermi to 
explain the acceleration of cosmic rays j^ Ulam formulated a mathematical 
model describing a massive particle bouncing between two infinitely heavy 
walls while one of the walls is oscillating [26] . A thorough analysis finally did 
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not fully confirm the expectations, however |2S1 [T21 [2U] . The model has also 
been reformulated in the framework of quantum mechanics [16] . 

More models of this sort have been studied in detail so far but we just 
mention one of them, the so called electron cyclotron resonance. One read- 
ily finds that electrons in a uniform magnetic field can gain energy from a 
microwave electric field whose frequency is equal to the electron cyclotron 
frequency. Because of an unlimited energy increase the relativistic effects 
cannot be neglected in a complete analysis. But even the relativistic model 
admits a quite explicit characterization of the resonant solution for a trans- 
verse circularly polarized electromagnetic wave propagating along the uni- 
form magnetic field [21j. In experimental arrangements the heated electrons 
are confined in a magnetic mirror field. Consequently, as they move along a 
flux tube of the mirror fleld they are exposed to the resonance heating only 
in a restricted region [2^ \TU\ [12] . This acceleration mechanism is widely used 
in plasma physics. 

Here we wish to discuss, on the quantum level, a model sharing some 
features with the preceding one. We again consider a charged particle placed 
in a uniform magnetic fleld. In our model the situation is simplifled, how- 
ever, in the sense that the particle is conflned to a plane perpendicular to 
the magnetic fleld. Instead of a transverse electromagnetic wave propagat- 
ing along the uniform fleld we apply, as an external force, an oscillating 
Aharonov-Bohm flux. The frequency of oscillations Q again coincides with 
the cyclotron frequency Uc or, more generally, it may be an integer multiple 

of Uc- 

The Aharonov-Bohm effect itself received a tremendous attention as a 
genuinely quantum phenomenon |2], and almost all its possible aspects have 
been studied in the time-independent case. For example, a careful analysis 
can be found in [22]. On the other hand, the time-dependent case represents 
an essentially more difficult mathematical problem and it has been treated 
so far only marginally in a few papers [TU [H [Sj |5] . 

The model we propose has already been studied in the framework of 
classical mechanics [1]. It turns out that a resonance acceleration again exists 
but it has some remarkable new features if compared to the standard electron 
cyclotron resonance. If Q is an integer multiple of Uc, then the classical 
trajectory eventually reaches an asymptotic domain where it resembles a 
spiral whose circles pass very closely to the singular flux line and, at the 
same time, their radii expand with the rate t^/^ as t approaches inflnity. The 
particle moves along the circles approximately with frequency Uc while its 
energy increases linearly with time. Denoting by £(t) the energy depending 
on time, an important characteristic of the dynamics is the acceleration rate 
which is computed in [4J and is given by the formula 

7acc:= lim^ = -^|$'(r)|. (1) 

t^oo t Air 

Here r is a real number which is expressible in terms of some asymptotic 



parameters of the trajectory. 

The purpose of the current paper is to demonstrate that one can derive 
a formula analogous to ([T]) also in the framework of quantum mechanics. To 
this end and because of complexity of the problem, we restrict ourselves to 
the case when the AB flux depends on time as a sinusoidal function. To 
this system we apply the quantum averaging method getting this way an 
approximate time evolution for which we observe a resonance effect whose 
principal characterization is again a linear increase of energy. 

Let us now be more specific. We consider a quantum point particle of 
mass M and charge e moving on the plane in the presence of a homoge- 
neous magnetic field of magnitude B. For definiteness, all constants M, e, 
B are supposed to be positive. Assume further that the particle is driven by 
an Aharonov-Bohm magnetic flux concentrated along a line intersecting the 
plane in the origin and whose strength $(t) is oscillating with frequency Vt. 

In the time-independent case, the Hamiltonian corresponding to a homo- 
geneous magnetic field and a constant Aharonov-Bohm flux of magnitude $o 
reads 

ff f 1 ^ ^ \ ( .^ e<l>o e5r2^^' 
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where (r, 6) are polar coordinates on the plane, and the Hilbert space in 
question is L^(]R_|_ x S"^, rdrd^). Making use of the rotational symmetry of the 
model we restrict ourselves to a fixed eigenspace of the angular momentum 
J3 = —ihdg with an eigenvalue j^h, J3 G Z. Put 

p:=j3-e%/{27ih). 

Then this restriction leads to the radial Hamiltonian 



Hip) = wi-A—drrdr + -,p + --] (2) 




in J^ = L^(]R_|_,rdr). Without loss of generality, we can assume that p > 
(note that H{—p) — H{p) is a constant). The boundary conditions at the 
origin are chosen to be the regular ones (then H{p) is the so called Friedrichs 
self-adjoint extension of the symmetric operator defined on compactly sup- 
ported smooth functions). Let us note that if < p < 1, then more general 
boundary conditions are admissible |7j but here we confine ourselves to the 
above standard choice. 
Let 

Uc = eB/M 

be the cyclotron frequency. The operator H[p) has a simple discrete spec- 
trum, the eigenvalues are 

En{p) = hujc{n + p + 1/2), n = 0,1,2,..., (3) 



with the corresponding normahzed eigenfunctions 
where 



2hJ \T{n + p+V 

are the normahzation constants and Ln are the generahzed Laguerre poly- 
nomials. 

Thus our main goal is to study the time evolution governed by the peri- 
odically time-dependent Hamiltonian H{a{t)) where 

a{t) =p + ef{nt) 

and f{t) is a 27r-periodic continuously differentiable function, fi > is a 
frequency and e is a small parameter. This means that the Aharonov-Bohm 
flux is supposed to depend on time as 

$(t) = $o-(27r;ie/e)/(fit). (5) 

Without loss of generality one can assume that 

f{t)dt = 0. (6) 

As discussed in [3J, for the values < p < 1 the domain of H{a{t)) in fact 
depends on t, and this feature makes the discussion from the mathematical 
point of view a bit more complicated. Nevertheless, the time evolution is still 
guaranteed to exist. 

2 The Floquet operator and the quasienergy 

Let U(t, to) be the propagator (evolution operator) associated with H{a(t)); it 
is known to exist [3] . An important characteristic of the dynamical properties 
of the system is the time evolution over a period which is described by the 
Floquet (monodromy) operator U{T,0), with T = 27t/Q. We are primarily 
interested in the asymptotic behavior of the mean value of energy 

{U{T,Of^,Hip)U{T,Ori;) 

for an initial condition ip as N tends to infinity while focusing on the resonant 
case when 

fl = fiUc for some /i G N. (7) 



A basic tool in the study of time-dependent quantum systems is the 
quasienergy operator 

K = -ihdt + H{a{t)) 

acting in the so called extended Hilbert space which is, in our case, 

jr = L^((0,T) xM+,rdMr). 

The time derivative is taken with the periodic boundary conditions. This 
approach, very similar to that usually applied in classical mechanics, makes 
it possible to pass from a time-dependent system to an autonomous one. The 
price to be paid for it is that one has to work with more complex operators 
on the extended Hilbert space. 

An important property of the quasienergy consists in its close relation- 
ship to the Floquet operator [TTl |27]. In more detail, if ^(t,r) G J^ is an 
eigenfunction or a generalized eigenfunction of K^ Kip = rjip, which also im- 
plies that ip{t + T,r) = ijj(t,r), then the wavefunction e~^^*^'^ip{t,r) solves 
the Schrodinger equation with the initial condition ipo{r) = il){0,r). It fol- 
lows that U{T, 0)'?/'o = e'^'^'^^'^'ipQ. Thus from the spectral decomposition of 
the quasienergy one can deduce the spectral decomposition of the Floquet 
operator. 

Let 

Ko = -ihdt + H{p) 

be the unperturbed quasienergy operator. Its complete set of normalized 
eigenfunctions is 

(here Z_|_ = {0, 1,2,.. .} stands for nonnegative integers, the wave functions 
iPn{p', r) are defined in (|1])), with the corresponding eigenvalues mhfl + En{p)- 
Thus Kq has a pure point spectrum which is in the resonant case ([7]) infinitely 
degenerated. 

To take into account these degeneracies we perform the following transfor- 
mation of indices. Denote by [x] and {x} the integer and the fractional part 
of a real number x, respectively, i.e. a; = [x] + {x}, [x] G Z and < {x} < 1. 
Furthermore, let 

be the remainder in division of an integer k by /i. The transformation of 
indices is a one-to-one map of Z x Z+ onto itself sending {m,n) to {k,i), 
with 

k = k{m, n) := fim + n, i = i{m, n) := [n/fi], (8) 

and, conversely, 

m = m{k, i) := [k/fj] - i, n = n{k, i) := fii + p(yU, k). (9) 



Using the new indices {k,i) we put 

^k,,{p;t,r)=T-'/'e''-^'''^'''Vn(k,i){p-^r). (10) 

Then the vectors "^k/, {k,i) G Z x Z_|_, form an orthonormal basis in the 
extended Hilbert space J^. For a fixed integer A; G Z let Pk be the orthogonal 
projection onto the subspace in ^ spanned by the vectors "^k/, ^ G Z+. Then 

Ko = ^XkPk where Xk = Hwdk + p + 1/2). (11) 



Furthermore, using the basis {'^k,e} one can identify J^ with the Hilbert 
space £^(Z X Z+). In particular, partial differential operators in the variables 
t and r like the quasienergy are identified in this way with matrix operators. 
In the sequel we denote matrix operators by bold uppercase letters. 

3 The quantum averaging method 

The full quasienergy operator K = K{e) depends on the small parameter e. 
Let us write K{e) as a formal power series, K{e) = Kq + eKi + e^K2 + . . .. 
In our case, 

K, = mt)fu.J^ + ^\^ K, = l^^, (12) 

■'^ ' \Mucr^ 2 J' 2Mr2 ' ^ ' 

and K^ = K4 = . . . = 0. The ultimate goal of the quantum averaging method 
in the case of resonances is a unitary transformation resulting in a partial 
(block- wise) diagonalization of K{e). Thus one seeks a skew-Hermitian op- 
erator W{e) so that e^^'^'>K{e)e~^'^'^^ commutes with Kq which is the same 
as saying that it commutes with all projections Pk- This goal is achiev- 
able in principle through an infinite recurrence which in practice should be 
interrupted at some step. Here we shall be content with the first order ap- 
proximation. 

Let us introduce the block-wise diagonal part of an operator A in J(f as 

dmgA:=^PkAPk. 

kei 

Thus diag A surely commutes with Kq. The off-diagonal part is then defined 
as offdiag A := A — diagA. Developing formally in e one has W{e) = eWi + 
0(e2) and 

Choosing Wi as 

Wi= Yl iK-K)-'Pk,K^Pk, 

k\,k2,ki^k2 



and 



one has 

[1^1,7^0] = - off diagi^i (13) 

and 

Let us note that the solution is also expressible in terms of averaging 
integrals, and this explains the name of the method [231 [IS]- In more detail, 
one has 

diagy4= lim - [ e-^"^°/''Ae^"^°/''du (14) 

Wi= lim -^ [ (r-M)e-'"^°/'^offdiag(iri)e'"^°/'"MM. (15) 

r->oo hr Jq 

After switching on the perturbation, any unperturbed eigenvalue A^ gives 
rise to a perturbed spectrum which, in the first order approximation, equals 
the spectrum of the operator XkPk + ^PkKiPk restricted to the subspace 
RanPfc C J(f . If the degeneracy of Afc is infinite then the character of 
the perturbed spectrum may be arbitrary, depending on the properties of 
PkKiPk. The corresponding perturbed (generalized) eigenvectors span a sub- 
space which is the range of the orthogonal projection 

Pfc(e) := e-^(^)p,e^(^) = P,-e[W^i,P,]+0(e2) 
= Pk-e{SkK,Pk + PuKA) + 0{^) 

where 

Sk = 2_^ ("^^ ~ '^^^ ^^ 

is the reduced resolvent of Kq taken at the isolated eigenvalue A^. Thus the 
first order averaging method is in fact nothing but the standard quantum per- 
turbation method in the first order but accomplished on the extended Hilbert 
space simultaneously for all eigenvalues of Kq (compare to [T71 Chp. II§2]). 
Our strategy in the remainder of the paper is based on replacing the true 
quasienergy i^(e) by its first order approximation 

ir(i):=i^o + ediagKi (16) 

and, consequently, U{T, 0) is replaced by an approximate Floquet operator 
f/(i) associated with -ft'(i). To determine the approximate Floquet operator 
[/(I) one has to solve the spectral problem for -ft'(i). To this end, as already 
pointed out above, one can employ the orthonormal basis {^E'fc^} in order to 
identify operators in J^ with infinite matrices indexed by Z x Z+. 

Let {e\;k G Z} denote the standard basis in £^(Z), and {e^;^ G Z+j 
denote the standard basis in £^(Z_|_). It is convenient to write ^^(Z x Z_|_) 
as the tensor product of Hilbert spaces £^(Z) (g) £^(Z_|_) which also means 



identification of the standard basis in P'iTL x Z_|_) witli tlie set of vectors 
{e\®e\\ keZJeZ+}. 

Let Pk be the orthogonal projection onto the one-dimensional subspace 
Ce^ C f{Z). Recalling ([IS]), dll]) and ([12]), the matrix ^(i) of the operator 
i^(i) expressed in the basis (ITP|) takes the form 

K(i) = ^Pfc®(Afc + eAfe) (17) 

fcez 

where Ak is the matrix operator in £^(Z_|_) with the entries 

{Ak)e,,e, = { ^kA,Ki^k,i2 )x ■ (18) 

To compute the matrix entries of Ak one observes that formally (see ([2])) 

K, = f{nt)dH{p)/dp (19) 

and so 

{^kA^Ki^k/^ )x = =^[/](^2 - ^i) {iPn{kM)ip)^ {dH{p)/dp) ^n(k/M) 
where 

p2it 

^[/](j) = (27r)-W e-'^-'fis)ds 
Jo 

stands for the jth Fourier coefficient of /. Recall that, by the assumption 

(Q, =^[/](0) = 0. Moreover, for ii ^ £2 one has n{k,ii) ^ ^(^,£2), hence 

{Vn{kA)iP)^H{p)ipn{k/2){p)) = 0- (20) 

In [3] it is derived that, for rii 7^ n2, 

{¥>nAP),^ ) = 777 r mm<^ — -,— -\ 21 

\ dp / 2{n2-ni) [liP^ni) 7(p;^2)J 

where 

'y{p;n) = [T{n + p+ l)/n\) 

Differentiating ( 120]) with respect to p and using ( |2T]) one finally obtains the 
relation 



^[/] (£2 - h) mini ^;^' "S^ tl? . ^i^' ''S^ !\\^ \ ■ (22) 



r4 1 =:^ 

^ "'''•'' 2^''^'^' "''-^^{^{p;n{k,h))'^{p;n{k,e,)) 

Note that n{k,i), as defined in ([9]), is yU-periodic in the integer variable 
k, and so is the matrix Ak, i.e. Ak+n = Ak- Moreover, since fiUc = 27r/T 
one also has e~^^''+i^^f^ = e"*-^*^/^* (see (ITT]) ). For an integer s, < s < ft, let 
J^s be the closed subspace in the original Hilbert space J^ = L^(M+,rdr) 



spanned by the vectors ips+jfi{i^), j = 0, 1, 2, . . .. Then J^ decomposes into 
the orthogonal sum 

jr = ^ © jTi © . . . © ^_i, 

and from the relationship between i^(i) and f/(i) , as recalled in Section |51 it 
follows that every subspace J^s is invariant with respect to [/(i). 

In the example which we study in more detail in the following section (for 
a sinusoidal function fit)), the matrix operators Ag have purely absolutely 
continuous spectra. For the sake of simplicity of the notation let us confine 
ourselves to this case. For a fixed index s, < s < /i, suppose that all gen- 
eralized eigenvectors and eigenvalues of Ag are parametrized by a parameter 
6 G (as,6s). Let us call them Xs{6) and i]s{6), respectively, i.e. 

AsXsie) = vs{e)xs{e), 

and write 

Xs{e) = iUid),Uiie),^sM^---)- 

The generalized eigenvectors Xs{6) are supposed to be normalized to the 6 
function, i.e. 

{x,{9,),Xsi92)) = 5i9i-92), 

which in fact means that ^s;e{9) as a function in the variables i G Z_|_ and 
9 G [as, bg) is a kernel of a unitary mapping between the Hilbert spaces 
£^(Z+) and L'^{{as,bs),d9). Thus the spectral decomposition of Ag reads: 

rbs 

yvefiZ+), Asv= / T]si9){x,{9),v)x,{9)d9. 

J a.. 



Put 



Then again. 



J2UA(^)^s+Up;r). (23) 



E.,{9ur)E,{92,r)rdr = 6(9,-92) 





and, for all ip^r) G J^, 

Uii)^{r) = e-2-(-+P+i/2)/A« / ^-^evs{e)T/H ^ s,(^), ^) S,(0, r) d^. (24) 



bs 



To get a correct approximation in the first order of the propagator one 
further has to take into account the transformation which is inverse to that 
generated by W{e) ~ e Wi. First observe that Wi is a multiplication operator 
on the Hilbert space J^ in the following sense. Let S be the unitary operator 
on J?^ acting as 

S^/j{t,r) = e'^'ij{t,r), W^P G JT. 
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An operator L on J^ commutes with S if and only if there exists a one- 
parameter T-periodic family of operators C{t) on L^(]R+,rdr) such that 
Li{j{t,r) = C{t)ip{t,r). Notice that 

S-^KoS = Kq + hVt. 

With this equality, it is obvious from (IT4l) that if A commutes with S then the 
same is true for diag A. Furthermore, as one can see from (TT2!) . Ki commutes 
with 5", and from fllSp one infers that Wi commutes with S as well. Hence 
there exists a one-parameter T-periodic family of skew-Hermitian operators 
Wiit) on L2(]R+,rdr) such that 

Wi%l){t,r) = >Vi(t)^(t,r), V^ G jr. 

Next notice that a transformation of the quasienergy operator of the 
form K = e^^*''i^e~^*^*\ where again W(t) is a T-periodic family of skew- 
Hermitian operators on L'^{M.^,rdr), implies a transformation of the associ- 
ated propagators according to the rule 

Hence the correct approximation of the Floquet operator reads 

UiT, 0) ^ f/app,o. = e-^^(°)t/(i)e^^^(°). (25) 

Let us note, however, that one has, for N &N and ip G J^, 

(f/,p%^,i7(p)[/,p^p,„,^) = (f/(i^^i,(if(p) + 6[>Vi(0),if(p)])^(i^^i)+O(e2) 

(26) 
where ipi = e'^^^'-'^V- If the commutator [Wi(0), H(p)] happens to be bounded 
then it does not contribute to the acceleration rate. 

Finally let us indicate how to compute the operator- valued function VVi(t). 
One has (here ipn = '^n{p',r)) 



oo 



^iW = E E e'''''Hj,ni,n,) (y.„„ ■ ) ^p^^ (27) 

j=—oo ni,n2=0 



1 -^ 



where 

W{j, m, ^2) = ^ / e"^^^* {ipn,,Wl{t)ipn,) dt. 

The commutator equation flT^ is equivalent to the differential equation 

-thWi'it) + [H{p),}Vi{t)] =offdiagiri. (28) 

Substituting ([27]) into (^ and using ([inD jointly with ([H]) one finds that 

^\f](j) I dH(n) \ 

w{j,ni,n2) = - — 7— ^\V^"i,^ V^na) a fij+ni-n2^0 (29) 

hUciiJ^j +ni- 712) \ op / 

and w{j,ni,n2) = otherwise. 
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4 A sinusoidally time-dependent AB flux 



In the remainder of the paper we discuss the example when f{t) = sin(t). 
The goal of the current section is to provide more details on the spectral 
decomposition of the averaged quasienergy i^(i) derived in (TTB]) . Naturally, 
rather than directly with the quasienergy we shall deal with its matrix, as 
given in f lT7|) and f lTSj) . 

We still assume that s G {0, 1, . . . , /x — 1} is fixed. For this choice of f(t), 
an immediate evaluation of formula fl2^ gives 



[-^s)ji,J2 



hjJr 



%2-ii|,isign(j2-Ji) n 



Atj< + s + iy 



1/2 



where j< = min{ji, J2}. Thus one has 

As = (hoojA) DJD^ 
where J is the Jacobi (tridiagonal) matrix with zero diagonal, 

\ 






ao 








ao 





ai 








ai 





a2 








a2 






v 



(30) 



and with the positive entries 



a. 



n 



yU j + S + Z/ 

\ M + S + P + 



1/2 



V 



and D is the unitary diagonal matrix with the diagonal (1, z, z^, i^, . . .). 

This is an elementary fact that the spectrum of J is simple since any 
eigenvector or generalized eigenvector is unambiguously determined by its 
first entry. Moreover, one readily observes that the matrices J and —J are 
unitarily equivalent, and so the spectrum of J is symmetric with respect to 
the origin. 

In our case. 



a, 



l-p/(2j) + 0(j-') as J ^00. 



Hence J is rather close to the "free" Jacobi matrix Jq for which aoj' = 1 for 
all j. The spectral problem for Jq is readily solvable explicitly (see below). 
It turns out that the spectral properties of J are close to those of Jq as well 
|13j . see also [22]. In particular, it is known that the singular continuous 
spectrum of J is empty, the essential spectrum coincides with the absolutely 
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continuous spectrum and equals the interval [—2, 2]. Furthermore, there are 
no embedded eigenvalues, i.e. if rj is an eigenvalue of J then |?7| > 2. 

Splitting J into the sum of the upper triangular and the lower triangular 
part, one notes that || J|| < 2 supjao, «i, a2, ■ ■ ■}■ In our example, aj < 1 for 
all j and so || J|| < 2 and, consequently, the spectrum of J is contained in the 
interval [—2,2]. This means that the only possible eigenvalues of J are ±2. 
But one can exclude even this possibility. In fact, suppose that Ju = 2u, 
with u = {uq, Ml, U2, ■ ■ .) and Uq = 1. Then 

ttj^iUj^i + ajUj^i = 2uj for j = 0, 1, 2, . . . 

(while putting u_i = 0). Summing this equality for j = 0,1, . . . ,n, and using 
that ttj < 1, one finds that Un+i > «„ + 1 forn = 0, 1, 2, . . .. Hence Uj > j + 1 
for all j, and so u is not square summable. Thus one can summarize that 
the spectrum of J is simple, purely absolutely continuous and equals [—2, 2 ]. 
Let us parametrize the spectrum of As = As{p) by a continuous param- 
eter 6^, < < vr, so that 



) 



ri{9) := {huJc/2) cos{9) 

is a point from the spectrum and x{p;6) is the corresponding normalized 
generalized eigenvector with components ^j{p] 0), j = 0, 1, 2, . . . (here we drop 
the index s at x and ^ in order to simplify the notation). The asymptotic 
behavior of the components C,j is known [TJl |6] ; one has 

e,(p; 9) ~ A{p; 6) V cos{j9 - (p/2) cot(^) log(j + 1) + 0(p; 9)) (31) 

for j ^ 0. Here A{p; 9) is a normalization constant and 0(p; 9) is a phase 
which depends on the initial conditions imposed on the sequence {C,j} (the 
initial condition is simply ^„i = 0) but the asymptotic methods employed 
in the cited articles do not provide an explicit value for it. In the limit case 
p = the generalized eigenvectors are known explicitly, namely 

^,(0;^) = V^^^sin((j + 1)^) 

for all J. Hence 0(0; 9) = 9 - 7i/2. 

The generalized eigenvectors are supposed to be normalized so that 

{x{p;9,),x{p;92)) = 6(9,-92). 

For p = 0, one can use the equality 



e^-- = nS{x) - V- 



oo 

1^' ■"^-' ■ 1-e- 

n=l 

which is valid for x = 9i — 92 & (~^7 ^) cind where the symbol V indicates 
the regular izat ion of a nonintegrable singularity in the sense of the principal 
value. The normalization is an immediate consequence of this identity. 
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For general p, the contribution to the S function should come from the 
most singular and, at the same time, the leading term in the asymptotic 
expansion of ^j{p;6), as given in fl3T|) . This time, when investigating the 
singularity near the diagonal 6i = 02 in the scalar product of two generalized 
eigenvectors, one is lead to considering the sum 



E 



laxinx 



where a = p/{2 sin^ ^i) is a real constant. Using the Lerch function <l>(z, s, v) 
one has for \z\ < 1 (see [91 § 9.55]), 



J^n'z"" = z<^{z,s,l) = T{1 - s) J2 (- log(^) + 27rm) 



-i+s 



n=l 



From here one deduces that, for any real a, 



E 



n'^^e''''' = 7t6(x) + iV- + gaix) (32) 

X 
n=\ 

where ga{p^) is a regular distribution, i.e. a locally integrable function. Hence 
in the general case, too, the normalization constant is given by 

As already mentioned, the phase 0(p; ff) in the asymptotic solution (13T]) 
remains undetermined. But we remark that a bit more can be said about the 
behavior of the phase near the spectral point (the center of the spectrum) 
which corresponds to the value of the parameter Q = 7r/2. More precisely, one 
can compute the derivative d(j){p; ti /2)/d6. Though this result is not directly 
used in the sequel it represents an additional information about generalized 
eigenfunctions of J . We briefly indicate basic steps of the computation in 
Appendix. 

5 The acceleration rate 

In the case when f{t) = sin(t) the commutator [>Vi(0), if (p)] occurring in 
fl2B]) can be shown to be bounded. This implies that instead of the approx- 
imate Floquet operator f/approx, as given in (^^, one can work directly with 
f/(i) defined in (12^ when deriving a formula for the acceleration rate. On 
the other hand, one should not forget about the transformation of the initial 
state, i.e. ^jJo has to be replaced by e^^^'^^^o, see (126|) . 

First let us shortly discuss the boundedness of the commutator. From (127]) 
and (12^ while using also (^1]) one derives that for ni,n2 G Z+, ni — ?22 7^ ±/i, 

4i {ni-Uif-iJ? XliP^ni)' 'j{p;n2) j' 
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Of course, the parallels to the diagonal determined by ni — n2 = ±/U can be 
explicitly evaluated as well but for our purposes it is sufficient to know that 
they are bounded. In [3], Lemma 6] it is shown that the matrix operator Q 
in £^(Z+) with the entries 

^"i>-2 4z (ni - n2) ""'"' 1 7(p; ^i) ' 7(p; ^2) J ^ ^ 

for ni 7^ n2 and otherwise is bounded. Thus to verify the boundedness of 
the commutator it suffices to show that the difference of matrices (J3S]) and 
fl34p has a finite operator norm. This can be readily done, for example, with 
the aid of the following estimate for the norm of a Hermitian matrix operator 
B [m § 1.4.3], 



\B\\ < sup ^ \Bn^ 



"iez+ „2=o 



"2 I 



To proceed further, we again fix an integer s, < s < /x. Suppose one is 
given a function q{0) G C^((0,7r)). Recalling (l23l) we put 



^(r)= / S,(^,r)f?(^)d^. (35) 

In what follows, we drop the index s and, whenever convenient, write simply 
H instead of H{p). Using ([2l]), one has, for A^ G N, 

^W, WTT^WA - f f J<cosei-cose2)uj,TN/2 



{U^i]^,HU^^)^) 



^0 



e 



x(S(^l),i7S(^2))f?(^l)^(^2)d^ld^2 



o-n •''0 



2 



j=0 

Note that {^j(p; ^)}j^o i^ ^^ orthonormal basis in L^((0, vr), d9) and so 

2 



E 







if^wrd^ 



Hence, in view of ([3]), the leading contribution to the acceleration rate comes 
from the expression 

fihu,f]{j + 1) r e-^"°«(^)-^™/2^^.(p; e)g{e) de 

Furthermore, restricting this sum to an arbitrarily large but finite number of 
summands results in an expression which is uniformly bounded in A^. This 
justifies replacement of ^j{p;9) by the leading asymptotic term, as given in 
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(PTj) (with A{p;6) = ^^/2/n). Hence the leading contribution to the acceler- 
ation rate is expressible as 

2hQ 



n 



^0 



h{9i, ^2)e^^('=°^^^-™^'^2)"'^™/2 ^(^^)^(^2) d^id^: 



2 



where 

oo 

h{e^, 62) = J2{j + 1) cos(j0i - I cot(0i) log(j + 1) + 0(p; e^ 

X cos(je2 - I cot(e2) log(j + 1) + (Pip; ^2)) . (36) 

The singular part of the distribution h{9i, 62) is supported on the diagonal 
61 = 92- The sum in fl36|) can be evaluated analogously as that in f l32|) with 
the result 

1 a_ 1 -YM^_iVw 



+ a regular distribution. 

Estimating the acceleration rate we can restrict ourselves to a sufficiently 
small but fixed neighborhood of the diagonal with a radius d > 0. Thus we 
arrive at the expression 

" U " U J- ^ -^ 

|6»i-6»2|<d 

Further we carry out the differentiation, as indicated in the integrand, and get 
rid of the terms which are not proportional to A^ or which are non-singular. 
Moreover, we use the substitution 6*2 = 9i + u. Thus we obtain the expression 



27r 



/ sm(6')|^^^^'^ 

TT Jo 



- sin [ - sm{e)ucTNu] dw J d^. 



Finally note that, for any a real, 

f^ 1 TT 

lim / — sinfaA^w) dw = — sign a. 

N^ooJq u ^ ' 2 ^ 



Suppose that the initial state is chosen as e~''^i'^°V- Then we conclude 
that the formula for the acceleration rate in the first-order approximation 
reads 

7acc := jim (f/5V^,i7(p)f/(^^>/(iVT||V^f) 

' ' ^ / sin(^)|^(^)|M0// W)\^dd. (37) 



/ ^0 
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Here we have used that 

nix 



Formula fl37p can be compared to formula ([T]), as derived for a classical 
particle, in the case when $(t) is given by ^ and f{t) = sin(t). Then ([T]) 
gives the acceleration rate 

7acc = \e\hwc^sm{^)/2 

where ^ G (0, tt) depends on some data which can be learned from the asymp- 
totic behavior of the classical trajectory. Let us finally note that, according 
to the analysis and discussion of the classical case presented in [4] , the first- 
order averaging approximation may in fact yield the correct acceleration rate 
(valid for the original system), and this is so even if the parameter e is not 
necessarily assumed to be very small. 



6 A numerical test 

We conclude our discussion by a presentation of a numerical result that 
agree quite nicely with the predicted acceleration rate (157|) . For the sake of 
simplicity we put Q = Uc = 1, and so /i = 1 and s = 0. We still assume that 
f{t) = sin(t). Concerning the physical constants, we set h = 1, e = 1 and 
M = 1. Furthermore, we choose p = 2.5, e = 0.4, and for the density g{6) 
determining an initial state according to fl35|) we take the Gaussian function 

/20\^/^ 
^(^)=(— J exp{-10{2-ef + 8te) 

restricted to the interval 6 G (0,7r). Its values near the limit points of the 
interval are in fact numerically indistinguishable from 0. Particularly, g{6) 
is normalized to unity with a negligible error, i.e. 

de = 1. 

The numerical method we use is based on expanding a solution of the 
time-dependent Schrodinger equation with respect to the time-dependent 
basis {(pn{a(t)); n G Z+}, with (pn{p) being defined in (j4]). Below we call the 
solution of the Schrodinger equation ipit)- Recalling fl23l) we put 

Mr)= [ Eo{e,r)g{9)d9, 
Jo 

and we have ||'?/^o|| = 1- The task is to solve the Cauchy problem for the 
time-dependent Schrodinger equation 

tdtm = H{a{t))tlj{t), ^(0) = tPo := e-^^^(°Vo. 
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Let us note that in the case of f{t) = sin(t) the matrix entries of Wi(0) are 
expressed as the finite sum 

(V5ni(p),>Vi(0)v3„2(p)) = w{l,ni,n2) + w(-l, rii, rig), 

with w{j,ni,n2) being given in ( 129|) (with /i = 1). To carry out the compu- 
tations we truncate the Fourier expansion of ip{t), 

oo 

^(t) = ^ Xnit)ipn{a{t)), Xn{t) = (v3„(a(t)), ^(t)), n = 0,l,..., 

n=0 

at some fixed order nmax- In this way we obtain a system of ordinary differ- 
ential equations for the Fourier coefficients 

^niax 

ix'^{t) = En{a{t))xn{t) -ia'{t) ^ {ipn{a{t)),ip'j{a{t))) Xj{t), 

j=0 

x„(0) = {ipn{a{0)),tjjo), n = 0, 1, . . . ,nmax- 

Explicit formulas for the scalar products are known from [3] (see (1211) ). In 
order to approximately solve this system we employ the explicit Runge-Kutta 
method of order 4 (RK4) with an adaptive step-size control, and we choose 

'^max iZU. 

From the computational point of view it is convenient to introduce the 
mean value of energy at time t as 

£it):={m,Hiait))m)- 
£{t) is then approximated by the sum 

£:(t)^5^i?„(a(t))|x„(t)|^ 

ra=0 

The acceleration rate is computed according to formula (13 7p in which one 
has to substitute ipQ for ip. Let us point out that this formula depends 
only on the time evolution over the intervals which are integer multiples 
of the period T, and clearly, H{a{NT)) = H{p) for A^ = 0,1,2,.... The 
predicted acceleration rate for the above particular values of parameters is 
7acc = 0.1796. The numerically computed function £{t)/t is compared to this 
value in Fig. 1. 

Appendix. The phase (^(p] 6) near the spectral 
point 

Here we compute the derivative 90(p; iT/2)/d9 of the phase 0(p; 9) introduced 
in fl3T]) . We know that always belongs to the spectrum of the Jacobi matrix 

17 



J introduced in (I5U1) . Putting u = {uo,Ui,U2, ■ ■ ■), with U2j+i = and 

u., = {-iy\i^^ (A.i) 

for j = 0,1,2,..., one has Ju = and Uq = 1. RecaUing that, in our 
example, aj = 1 — p/{2j) + 0(j^^) one derives that 

U2j = (-l)%oo(l +p/(8j) + 0(j-2)) as J -^ oo, 

where 

Moo = lim {-iyu2j 

is a finite constant (depending on p, however). Comparing to flM]) . with 
A(p; 0) = A/2/vr and = 7r/2, one finds that 

x{p; 7r/2) = {^/2/tt/uoo) u. 

Moreover, 0(p; 7r/2) = 0. 

Differentiating the equahty 

Jx{p;9) = 2cos{9)x{p;9) 

with respect to 6 at the point 7r/2 and using the substitution 

dx{p;7i/2)/d9 = - (2^/2/^ /ujjv, 

with V = (t>o, f 1, ^2, • • •); 0^6 arrives at the equation Jv = u. From (13T|) one 
deduces that 

,, ~ i .„ .„(, ^) (, + I logo + 1) + '-^^) (A.2) 

for j ^ 0. This suggests that one can seek a solution v such that V2j = 
for all j. This assumption on v is in fact necessary and makes the solution 
unambiguous since otherwise one could add to v any nonzero multiple of u 
thus violating the asymptotic behavior (IA.2p . Given that all odd elements 
of the vector u and all even elements of v vanish the equation Jv = u 
effectively reduces to a linear system with a lower triangular matrix which is 
explicitly solvable. Using flA.ip one can express the solution as 

1 ^ 

^2,+i = — y](«2fc)', J =0,1,2,.... (A.3) 



''■ "2,^2, ,^, 



Noting that 






fc=0 
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where 7^; is the Euler constant, and that 

{u,,f = ui{i+p/m+o{k-')) 

one derives 

,_„ \"oo/ "± 



fc=0 



00 / / \ 2 



+ Em^ -i-s^ -^oir-; 



fc=0 \ ^ °°/ 



Using (1A.3P and comparing to (1A.2P one finally arrives at the relation 



l + |(M2)+,.)+2g((g^l 



•9^ 2^ "^ ^ '"^ ^VVwoo; 4j + 2 
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Figure 1. The function £{t)/t, with £(t) = {ijj{t), H{a{t))^{t)) and ^(t) being 
a normahzed solution of the time-dependent Schrodinger equation, compared 
to the value of the acceleration rate 7acc derived in eq. flTTI) . 
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